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Precipitation in aqueous mixtures with addition of strongly hydrophilic or 
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We examine phase separation in aqueous mixtures due to preferential solvation with a low-density 
solute (hydrophilic ions or hydrophobic particles). For hydrophilic ions, preferential solvation can 
stabilize water domains enriched with ions. This precipitation occurs in wide ranges of the tem- 
perature and the average composition above a critical solute density n p , where the mixture solvent 
would be in a one-phase state without solute. The volume fraction of precipitated domains tends to 
zero as the average solute density n is decreased to n p or as the interaction parameter y is decreased 
to a critical value Xv K we start with one-phase states with n > n p or y > y p , precipitation pro- 
ceeds via homogeneous nucleation or via heterogeneous nucleation, for example, around suspended 
colloids. In the latter case, colliod particles are wrapped by thick wetting layers. We also predict a 
first-order prewetting transition for n or y slightly below n p or y p . 

PACS numbers: 82.45.Gj, 61.20.Qg, 64.75. Cd, 81.16.Dn 



I. INTRODUCTION 

In fluid mixtures composed of water and a less polar 
organic liquid, phase separation behavior can be dras- 
tically changed by a small amount of salt QUI]- With 
addition of a 10~ 3 mole fraction of hydrophilic salt like 
NaCl, the coexistence curve is typically shifted by a few 
Kelvins upward in the USCT case or downward in the 
LSCT case. In some aqueous mixtures, even if they are 
miscible at all T at atmosphere pressure without salt, 
addition of a small amount of salt gives rise to reentrant 
phase separation behavior [3j-|6| . 

With ions added, moreover, many groups have also 
found long-lived heterogeneities (sometimes extending 
over a few micrometers) in one-phase states [9( and a 
third phase visible as a thin solid-like plate at a liquid- 
liquid interface in two-phase states [l(| • A representative 
system is a mixture of H2O (or D2O), tri-methylpyridine 
(3MP), and NaBr. The heterogeneities have been de- 
tected by dynamic light scattering [111 ], indicating spon- 
taneous formation of aggregates or domains with small 
diffusion constants. Their size and the volume fraction 
arc 50A and 10" 6 as a typical example. The third phase 
have been observed also in isobutylic acid (IBA)-water 
mixtures without added salt, where the isobutyric acid 
partly dissociates into Butyrate - and H + ions. In ad- 
dition, a drastic decrease in the ion mobility was ob- 
served in IB A- water for high contents of IBA [l^. Re- 
cently, mesophases have been observed when sodium 
tetraphenylborate NaBPli4 was added to D2O-3MP at 
about 0.1 M [l3[ . This salt dissociates into hydrophilic 
Na + and hydrophobic BPI14 . The latter ion consists of 
four phenyl rings bonded to an ionized boron. 

Dramatic ion effects are ubiquitous in various soft ma- 
terials such as polymers, gels, colloids, and mixtures con- 
taining ionic surfactants. For example, complex effects 
are known to be induced in polyelectrolytes when a sec- 
ond fluid component (cosolvent) is added to water. In 
particular, precipitation of DNA has been widely ob- 



served with addition of alcohol such as ethanol to water 
14{. Here the alcohol added is excluded from condensed 



DNA, which suggests solvation-induced wetting of DNA 
by water. We also mention observations of crystal for- 
mation of micron-sized water-rich droplets in a less polar 
oil-rich phase with addition of HBr [15| . The mechanism 
of this ordering was ascribed to the Coulomb interaction 
due to asymmetric partitioning of cations and anions in 
water and oil [ill ]. 

In these phenomena, the solvation (or hydration) in- 
teraction among ions and polar molecules should play a 
major role. Sma ll hy drophilic ions are solvated by several 
water molecules [17H19| and the resultant solvation chem- 
ical potential per ion much exceeds the thermal energy 
ksT and strongly depends on the ambient composition 
of the solvent. On the contrary, some hydrophobic parti- 
cles dislike to be in contact with water molecules, which 
tend to form aggregates in water and are more soluble 
in oil than in water. In aqueous mixtures, a hydrophobic 
ion should be surrounded by oil molecules, which explains 
the mesophase formation observed by Sadakane et al [l3| ■ 
Hydrophilic and hydrophobic particles strongly affect the 
surrounding hydrogen bonding network in different man- 
ners [17( and the solvation is highly cooperative. Fur- 
thermore, when a small amount of water was added to 
methanol-cyclohexane, the coexistence curve was largely 
shifted [20| and a water-induced, methanol-rich wetting 
layer appeared (which was nonexistent without water) 
[2l| . These findings indicate strongly selective molecular 
interactions between water and the two components. 

Obviously, the preferential solvation should strongly 
affect the phase transition behavior in aqueous mixtures 
with a small amount of hydrophilic or hydrophobic so- 
lute, though this aspect has not been well studied. Re- 
cently, some efforts have been made to elucidate the sol- 
vation effects in phase separation in mixture solvents in 
electrolytes (22b 28] . polyelectrolytes [29j . and ionic sur- 
factant s 1301] . A review on the static properties was pre- 
sented [31]. In the dynamics, a number of problems still 
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remain unexplored [271 12SJ ■ 

In this paper, we present a theory of solvation-induced 
phase separation with hydrophilic or hydrophobic so- 
lute. Thus we aim to explain the observed heterogeneities 
in aqueous mixtures [9J. We may suppose hydrophilic 
monovalent ion pairs such as Na + and CI - in a binary 
mixture of water and a less polar component (called 
oil). As is well-known, hydrophilic ions induce clus- 
tering of water molecules around them on microscopic 
scales, forming a solvation shell 

EMU- 

We shall see 

that strongly hydrophilic ions can moreover induce for- 
mation of water-rich domains on macroscopic scales for 
sufficiently strong preference of water over oil. This can 
occur even when the mixture is outside the coexistence 
curve without ions. 

The organization of this paper is as follows. In Sec. II, 
two-phase coexistence induced by the preferential solva- 
tion will be studied numerically and theoretically, where 
the electrostatic interaction does not appear explicitly. In 
Sec. Ill, inhomogeneous profiles such as interfaces will be 
calculated. We will show that a precipitated droplet can 
be stable only above a minimum radius due to the surface 
tension effect and that a solvation-induced prewetting 
transition occurs far from the coexistence curve without 
solute. In Sec. IV, two-phase coexistence with hydrophilic 
ions will be examined in the presence of the charge effect. 
In Sec.V, we will investigate the precipitation from one- 
phase states taking place as homogeneous nucleation. In 
Sec. IV, we will investigate adsorption and precipitation 
on colloid surfaces. 



II. PHASE SEPARATION WITH STRONGLY 
SELECTIVE SOLUTE 

Neglecting electrostatic interaction but accounting for 
the solvation interaction, we first consider a binary mix- 
ture composed of water and a less polar component in a 
cell with a fixed volume V. The second component will 
be simply called oil hereafter. Ions will be treated as a di- 
lute third component (solute) with density n(r). We may 
also suppose strongly hydrophilic or hydrophobic neutral 
particles (possibly with a complex structure) added in an 
aqueous binary mixture. Such particles can be solvated 
by a certain number of water or oil molecules depending 
on whether they are hydrophilic or hydrophobic. 

In the following, we will suppose hydrophilic ions or 
particles. However, our results can be equally applica- 
ble to hydrophobic ions or particles if water and oil are 
exchanged (or <p is replaced by 1 — 0). The Boltzmann 
constant fcs will be set equal to unity. 



A. Conditions of two phase coexistence with solute 

The volume fractions of water, oil, and solute are writ- 
ten as 4>(r) and (j>'(r), and Ujn(r), respectively, where vi 
is the solute volume. If the two solvent species have the 



same molecular volume vo, their densities are (j>/vo and 

1 /3 

4> /vq. The solvent diameter is of order a = v . The 
space-filling condition is written as 

(j> + 4>' + vin = 1. (2.1) 

The solute volume fraction is assumed to be small or 
V]n <C 1. This is more easily satisfied if the solute size is 
smaller than that of solvent or vj < vq. In this paper, to 
simplify the calculations, we thus set 

<t)' = 1 - cj). (2.2) 

The free energy density /tot {4>, n) for the composition 
<f> and the solute density n consists of three parts as 

/tot (6 n) = f(4>) + TnHnX 3 ) - Tgntp. (2.3) 

The first term is assumed to be of the Bragg- Williams 
form for a fluid mixture [32j , 

f(4>) = - + (1 - 0) ln(l -4>) + - 4>)], (2-4) 

vo 

where x is t ne interaction parameter dependent on T. 
1/3 

The length a — v a represents the molecular diameter 
(~ 3A for water). In Eq.(2.3), the second term arises 
from the solute translational entropy. The A is the ther- 
mal de Broglie length (but the terms linear in n are irrele- 
vant and A 3 may be replaced by vq in the following). The 
third term arises from the solute preference of water over 
oil. The parameter g is assumed to much exceed unity 
(;§> 1) 23]. We fix the average densities of the constituent 
components, which are expressed as 

4> = J dr<j)/V, n = J drn/V. (2.5) 

In our theory, Xi 4>i an d n are relevant control param- 
eters. Changing \ through the critical value \ c = 2 is 
equivalent to changing T through the solvent critical tem- 
perature T c . 

In two phase coexistence in equilibrium, let the com- 
position and the solute density be (4> a ,n a ) and (4>p,np) 
in phases a and /3 with (ftp < <f> < (j) a and np < n < n a . 
Here we give simple thermodynamic arguments. We 
define the chemical potentials as h — dftot/dcf) and 
M = dftot/dn. From Eq.(2.3) we obtain 

h = /'(0) - Tgn, (2.6) 
,i = T[ln(nA 3 ) + l- 5 0], (2.7) 

where /' = df/<f>. First, the homogeneity relation 
//(0 a , n a ) = fi((pp,np) for solute gives rise to 

n a = A e 9 ^, np = A e 9 ^. (2.8) 

The coefficient Aq = A~ 3 exp(/z/T— 1) is determined from 
the conservation of the solute number in Eq.(2.5) as 

A = fi/{e?*), (2.9) 
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where (■ • •) denotes the space average in the cell. Let 7q 
be the volume fraction of phase a. Neglecting the volume 
of the interface region, we obtain 

(e»*) = 7ct e 9 ^ + (1 - la )e 9 ^ (2.10) 

From Eq.(2.5) "f a is expressed in terms of 4> & n d n as 

= (n- rip) /An. (2.12) 

The differences between the two phases are written as 

A0 = (j) a - (j>p > 0, An = n a - np > 0. (2.13) 

With the aid of Eq.(2.8), the ratio n a /n is written in 
terms of the compositions as 

n a /n = A(/>/[$ ~ fa) + (&» ~ ^)e' a ^]. (2.14) 

For not small A(j>, we find n a /np — e gA ^ ^> 1 and An = 
n a . Namely, the solute density is much higher in phase a 
than in phase (3. The solvation part —Tgn a 4> of the free 
energy density is significant in phase a even for small h. 

Next, the homogeneity relation h(<fi a ,n a ) = h{<pp,np) 
for the solvent composition is written as 

h = f'(<f> a ) ~ Tgn a 

= f'(<j>p)-Tgnp. (2.15) 

In equilibrium, we also require minimization of the grand 
potential density u> defined by 

w = /tot - h(f> — fin 

= f-hcj)-Tn, (2.16) 

where the second line follows from Eq.(2.7). In the two 
phases, u assumes the same value so that 

f{<fia) ~ fi't'p) ~ TAn = hA(j). (2.17) 

We may derive Eqs.(2.15) and (2.17) from minimiza- 
tion of the total free energy F under Eq.(2.5). Using 7 a 
we express F as 

F/V = J a ftot(<Pa,n a ) + (1 - la)ftat(<l>p,np), (2.18) 

where the surface free energy is neglected. If use is made 
of the relation n(ln(nA 3 ) — g<p) — nlnAo in the two 
phases, the above expression is rewritten as 

F/V = 7a/(0 Q ) + (l-7a)/(^) + Tnln(nA 3 ) 

- Tnln[ 7ct e^° + (1 - 7 Q )e 9 ^]. (2.19) 

The third term in the first line is a constant at constant 
n = (n) and is irrelevant, but the fourth term in the 
second line is a singular solvation contribution at fixed 
n = (n) relevant for large g (even for small n). Since 
<f> = (4>) is also fixed, we should minimize 
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FIG. 1: (Color online) Left:Compositions 4> a and <f>@ vs %■ 
Right:Semi-logarithmic plots of volume fraction 7^ and nor- 
malized ion density n a /n of the water-rich phase a vs \. 
Here n = 6 x lO -4 ^ 1 and g — 11. Precipitation occurs 
for Xp < X i; 2. For <j> = 0.35 (top), <j> a changes continuously. 
For <j) = 0.46 (bottom), (f> a jumps at x — 2. Plotted in the 
left also are the coexisting region without solute (in blue) and 
the spinodal curve with solute (see Eqs.(3.7) and (5.2)) [23| . 



with respect to </> Q , (ftp, and 7ct . Here h appears as 
the Lagrange multiplier. With the aid of the expression 
(2.19), the minimum conditions, dui/d<p a — duj/dtpp = 
dCj/dj a = 0, readily lead to Eqs.(2.15) and (2.17). 

Without solute n = 0, two-phase coexistence is pos- 
sible only for x > 2. Thus, in the appendix A, we will 
perform the Taylor expansions of <j) a and with respect 
to n for x > 2. We shall see that the water-rich coex- 
istence branch is much deformed even by a very small 
amount of highly preferential solute (as can be seen in 
Fig.l below). 



B. Numerical results of precipitation 

Next we give numerical results on the phase behavior 
of our system. We will set g = 11 mostly, but will also 
present additional results for g — 15 in Fig. 3 below. The 
solute density n will be measured in units of Vq 1 . 

In Fig.l, we display <j) a and <f>p in the left panels and 
n a and j a in the right panels as functions of x ( or T 
experimentally) at n = 6x lO -4 ^ 1 . Due to the nonlinear 
solute effect, a precipitation branch appears in the range, 



u> = F/V - hfr a (p a + (1 - 7a )^ - $\ (2.20) 



X P (0,n) < x < 2, 



(2.21) 
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FIG. 2: (Color online) Volume fraction 7^ (top) and com- 
position 4> a (bottom) of the water-rich phase a in the x~4> 
plane at fixed n (left) and in the n-(f> plane at fixed \ (right) 1 
where g = 11. Precipitation occurs for \ > Xp(0>^) (left) 
and n > n p (4>, \) (right) outside the slashed regions. In the 
left bottom, a first-order transition occurs for \ = 1.95 and 
4> > 0.395 (in red). See the bottom plates in Fig.l. 



where the solvent would be in one-phase states without 
solute. Here 7 Q tends to zero as x ~ ► Xp> where the lower 
bound Xp depends on <fi and h. This branch appears un- 
der the condition e 9 ^'^ 1 (see the next subsection). 
For large g{= 11 here), the precipitated domains become 
solute-rich (salted water for hydrophilic ions) with 



b a = 1, n a > np — e 9 ^n c 



(2.22) 



The right panels are on semi-logarithmic scales. With 
decreasing x, we can see an increase in the solute density 
n a and an decrease in the volume fraction j a . Details 
of the figure are as follows, (i) For (f> = 0.35 (top), (f> a 
changes rather abruptly around x ~ 2 but continuously 
as a function of x, being minimum at x — 2.05. The 
n a increases up to D.381^ 1 as \ — > Xp = 0.687. (ii) 
For <fi ~ 0.46 (bottom), the precipitation branch much 
shrinks with Xp = 1-834, while n a — 0.0842w ( ^ 1 at x = 
X P . For this larger 0, 4> a changes discontinuously on a 
hysteresis loop in the range 1.953 < x < 1-967. If F in 
Eq.(2.18) is minimized, a first-order transition is found 
to occur at x = 1-957. Even for x > 2, the two-phase 
coexistence is much deformed by the solute. 

In Fig. 2, we show 7 Q and <pa hi the x~4> plane at h = 
6 x lO -4 ^ 1 (left) and in the h-(f> plane at x = 1-8 (right). 
For fixed x m the n-0 plane, precipitation occurs for 

n>n p {4>, X ), (2-23) 

where 7^ tends to zero as n — > n p . This minimum solute 
density n p depends on 4> and x- As the crossover to the 
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FIG. 3: (Color online) Left:x P (0, n) vs <j> for three values of n 
at g = 11 (top) and 15 (bottom), which nearly coincide with 
the asymptotic formula (2.33) (dotted line) for <f> < 0.35 and 
converge to the coexistence curve for larger <j>. Coexistence 
region without ions is also shown (in blue). Right: von p ((j), \) 
vs on a semi- logarithmic scale for three values of x at g = 
11 (top) and 15 (bottom), It coincides with the asymptotic 
formula (2.34) (dotted line) for x < 2 and tends to zero at 
the coexistence composition, 0.249 (top) or 0.204 (bottom). 



asymptotic behavior in Eq.(2.22), <p a becomes apprecia- 
bly smaller than 1 for x ~ 2 or for cf> ~ 0.5. For large 
cj) > 0.5, the precipitation branch shrinks to vanish for 
small solute densities, where water molecules are already 
abundant around solute particles in one-phase states. 

In Fig. 3, at g = 11 and 15, we show Xp(^ ; n) vs 4> f° r 
three values of n (left) and n v {4>, x) v s 4> f° r three values 
of x (right). These quantities decrease dramatically at 
small <p and their magnitudes strongly depend on g. Even 
for x > 2, n p can be calculated outside the water-rich 
branch of the coexistence curve and it tends to zero as <j> 
approaches the coexistence composition without solute. 
Theoretical explanations of their behavior will be given 
in the following subsections. 



C. Theory of asymptotic behavior 

We present a theory on the asymptotic behavior of the 
precipitation branch for g 3> 1 in the region x < 2. At its 
starting point, we assume the branch satisfying Eq.(2.22) 
and confirm its existence self-consistently. Since (f> a = 1, 
the logarithmic term (oc (1 — </>)ln(l — 4>)) in the free 
energy density (2.4) is crucial in phase a. 

We first neglect the term —Tgnp in Eq.(2.15) from 
gvonp <C 1 and the term f(<p a ) in Eq.(2.17) from f(<fi a ) ~ 
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FIG. 4: (Color online) v G(cj>)/T denned in Eq.(2.26) in the 
4>~X plane. It is negative only in the inner green region, where 
the minimum of \ is 2.455 at tj> = 0.285. It is positive outside 
it in the blue region. The coexistence curve without solute is 
written in red, on which the minimum of x is 2 at = 0.5. 



T(l - 4> a )\n(l - <j) a ). Note that gvnnp ~ 1(T 2 in Fig.l. 
It follows a simplified equation, 

h = f'ifa) = -[/(^) + Tn a ]/(1 - M. (2.24) 

This determines the solute density n a in phase a as a 
function of <f)p in the form, 

n a = G(4>p)/T, (2.25) 

where G(4>) is defined by 

G(0) = -ftf) - (l - 4,)f{ff>) 

= -(T/«o)M + x(W) 2 ]- (2-26) 

Using the free energy (2.4), we obtain the second line and 
plot it in the <j>-x plane in Fig. 4. It is positive outside the 
coexistence curve, ensuring n a > in Eq.(2.25). 

In Eq.(2.15), we next use Eq.(2.4) for f'((f> a ) to obtain 

h £* »„ *T[- ln(l - <f> a ) - X ] ~ Tgn a = f'fo), (2.27) 

where the logarithmic term (oc ln(l — cj>)) balances with 
the solvation term (oc gn a ). Use of Eq.(2.26) gives 

l — </>n = Ap exp[-gvoG(if>p)/T\, (2.28) 

where the coefficient Ap is given by 

A^=exp[-x-uo/'(^)/T], (2.29) 

so A/3 is of order unity. The factor exp[—gvoG(<p/3)/T] in 
Eq.(2.28) is very small for j > 1, leading to <f) a = 1 as 
in Figs.l and 2. 

Furthermore, from Eqs. (2.12) and (2.25), the volume 
fraction "/ a of phase a is approximated as 

la = n/n a - e~ 9 ^ = Tn/G(4>p) - e~ 9 ^ . (2.30) 

The above relation is rewritten as 

Gfo) s Tn/{ 7a + e' 9 ^) 



4> - 0/3 + (1 - <Pp) exp[-g(l - <j)p) 



From the first to second line, we have used Eq.(2.11) 
and replaced A</> by 1 — c/>p. This equation also follows 
from Eqs. (2. 14) and (2.25). It determine (ftp and -f a = 
{4> — (f)p)/{l — (j)p). We recognize that G{<j)p) increases up 
to Tne 9 ^ Tne 9 ^-^ as j a -> or as <f>p -t 4>. In this 
limit it follows the marginal relation, 

G{$) = Tne 9 ^-® (j a -> 0). (2.32) 

If n is fixed above n p , Eq.(2.32) holds at x — Xp so that 

Xp = [-\n4>-v Q ne 9 ^}/(l~4>) 2 , (2.33) 

where we use the second line of Eq. (2.26). Notice that the 
solute density n appears in the combination Se s ' 1_< ^ and 
its effect of lowering Xp is much amplified for g{l — <j)) 3> 1 
even for very small n. (For n — the right hand side of 
Eq.(2.33) is above the coexistence curve.) On the other 
hand, if x is fixed above x P , Eq.(2.32) holds at n = n p . 
Thus the minimum solute density n p is estimated as 

n p = e~- 9{1 -^G(0)/T, (2.34) 

which is much decreased by the factor e -9 ^ 1 ^ . 

In Fig. 3, the curves of % p and n p nearly coincide with 
the asymptotic formulas (2.33) and (2.34) in the range 
<f> < 0.35 for Xp an d m the range x < 2 for n p . They 
exhibit a minimum at ip ~ e~ 9 /v^ng for Xp an d at 4> ~ 
g _1 for n p . For larger (j> > 0.35, % p nearly coincide with 
the coexistence curve, indicating disappearance of the 
precipitation branch. Notice that n p decreases to zero 
as <j> approaches the coexistence composition (j> cx = 0.249 
at x = 2-2 (top) and 0.204 at % = 2 -3 (bottom), where 
phase separation occurs without solute. 

D. Theory for small volume fraction of precipitates 

We next construct an analytic theory of precipitation 
for small volume fraction 7 a « 1. We will calculate F 
in Eq.(2.19) treating j a as an order parameter of the 
phase transition. Though this theory is applicable only 
for 7 Q <C 1, it can yield Eq.(2.21) in the asymptotic limit. 

In Eq.(2.19) we first expand (j>p and f((f>p) with respect 
to 7 a as 

<f>0 = (<!>- 0a)/(l - la) = 4> + (4> - <t>aha, 
/(^)=/(0)+/'(0)(0-^a)7a. (2.35) 

Let F = Vftot{(t>,n) be the one-phase value of F (at 
7a = 0). The increment AF = F — Fq is due to precipi- 
tate formation and is of the form, 

AF/V = n((t> a ) la - Tnln[l + *(^) 7q ], (2.36) 

where uJb(4>) an d ^(4>) are defined by 

tt(0) = f{4>) - - - $), (2.37) 

= exp[ 5 (0 -$)}-!- g(<j> - $), (2.38) 
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FIG. 5: (Color online) fi(0)/*(0) multiplied by 10 4 u /T for 
three values of \ at g = 11 and = 0.46. Each curve exhibits 
a minimum at = (j> a , where <j> a is estimated as in Eq.(2.51). 
Each minimum value is equal to 10 4 i>on p from Eq.(2.46). 
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FIG. 6: (Color online) Theoretical <f> a in the range [0.975, 1] 
(left) and n p multiplied by 10 3 «o (right) for <j) = 0.25, 0.35, 
and 0.45 at g = 11. These curves are calculated for j a <C 1 
from Eqs.(2.45) and (2.46). For these parameters, they only 
slightly deviate from the asymptotic formulas in Eqs.(2.49) 
and (2.52). The curves of n p are nearly linear in \ from the 
linear dependence of G(<j>) on \- 



where (f> is treated as a constant and its dependence 
of £1 and $ are suppressed. In Eq.(2.36) we have not 
expanded the logarithmic term with respect to j a , be- 
cause the coefficient ^>(4> a ) grows strongly as e 9 ^ a ~^ 
for g(4> a — 4>) 3> I- For x < 2, the positivity fl((f) a ) > 
follows in the region (f> < <j) < <t>a since Eq.(2.37) gives 



(2.39) 



where /"(<£) = d 2 f/d 2 <j) > 0. Note that if 4> and cj) in 
£l((f>) are replaced by 1 and </>p, respectively, fi(</>) be- 
comes identical to G((t>p) in Eq.(2.26). If use is made of 
Eq.(2.4), ft(cf>) is explicitly written as 



O-x(0-0) 2 - (2.40) 



^n(0) = 0in(|)+(i-0)Mf 



For each given 4> and n, AF should be minimized with 
respect to 7 Q and <j) a . Here (ftp is given by Eq.(2.35). 



From Eq.(2.36) its derivative with respect to 7 Q is 

1 dAF Tm{<t> a ) 

The right hand side can vanish only when 

Z(cj> a ) ee fi(0 a )/Tn*(0„) < 1. (2.42) 

If Z{(j) a ) < 1, j a is nonvanishing and is written as 

la = Tn/Q(4> a ) - (2.43) 

For this j a , AF in Eq.(2.36) becomes negative as 

AF/V = Tn[l- Z((f> a )+hxZ(4> a )]. (2.44) 

In the range < Z < 1, the function 1 — Z + \nZ is 
negative and decreases with decreasing Z. In order to 
minimize AF, we should thus minimize Z((f>) or the func- 
tion f2(0)/^(0) at 4> = <f) a . The equation to determine 
4>a = 4>a{<i>, x) is hence given by 



d n(0 a ) _ Q'(0„) fi(^ Q ) 



90 Q (0 Q ) *(0 a ) *(^ a ): 



= 0, (2.45) 



where Q'(cf>) = dVt(4>) / d4> and *'(^) = d$((f>)/d(f>. Note 
that (/> Q is independent of n for "f a <C 1. In Fig. 5, we plot 
fi (</>)/ "J (</>) for three values of \ a * = 0.46 to demon- 
strate the existence of its minimum at </> = <j) a . In the 
left panel of Fig. 6, we show Q calculated from the above 
equation. 

In terms of the above <p a (dependent on (j) and x), the 
minimum solute density n p is obtained from the condi- 
tion Z((j) a ) = 1 as 



n p = fi(0 Q )/r*(0 Q ), 



(2.46) 



which is a function of <p an d X- Then Z((j) a ) — n p /n. 
The condition Z((f> a ) < 1 yields n > n p , as ought to be 
the case, and ^ a is also expressed as 



7« = (n/n p - 1)/*(0 Q ). 



(2.47) 



In the right panel of Fig. 6, we show n p calculated from 
Eq.(2.46) using the data of <fi a in the left panel. 

After <p a and n p have been determined, Eq.(2.47) indi- 
cates that an arbitrary value can be assigned to n in the 
range < n/n p — 1 -C ^((fia). Furthermore, by decreas- 
ing x &t fixed n and <f>, we can make Z(4> a ) approach 
unity, where x ~ * Xp = Xp(^j^) m Eq.(2.21). Using 
Eq.(2.40) we may express Xp as 



<P 



111 



1 



1 - 



(2.4 



The relations (2.46) and (2.48) for the lower bounds are 
exact in our model and are consistent with the asymp- 
totic ones (2.33) and (2.34) in the limit of (f) a = 1 and 
9(A a ) = e«0--4>) > 1. 
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We now give analysis of the function il(<j>) / ^ ((/)) in the 
range <fi < <f>. If ^ is not very close to unity, f2(</>) /^(4>) de- 
creases rapidly with increasing (j>, because l/^>((f>) oc e~ ff ^ 
decreases rapidly. However, as 4> — > 1, the logarithmic 
term in f(<p) eventually comes into play. To show this, 
we approximate £l(<j>) for <j> = 1 as 



n($<*G($ + -(i 

«0 



In 



1 



(2.49) 



where G{<j>) is defined in Eq.(2.26) and ln(A ) = — x — 
v f'{<j>)/T_. The coefficient Ap in Eq.(2.29) tends to A 
as c/)/3 — > (j). The derivative of fi((/>)/\&(0) with respect to 
6 is thus estimated as 



d~4>M^) 



— In— — 

v A 



-g(<t>-4>) 



(2.50) 



In this derivative, the logarithmic term grows weakly but 
can balance the solvation term (oc g) as <fi — > 1, so that 



1 — <§>a — Aq CX P (- 5UO G(0)/T) 



(2.51) 



in accord with Eq.(2.28). The asymptotic form of n p 
is written as in Eq.(2.34) if we set Q(4> a ) — G(<f>) and 
^ e» A * in Eq.(2.46). 



III. INHOMOGENEOUS COMPOSITION 
PROFILES 

Here, still neglecting the charge effect, we examine the 
surface tension, the stability of a spherical droplet, and 
the surface adsorption near a boundary wall. These prob- 
lems stem from the solvation-induced phase separation. 



A. Gradient free energy and surface tension 

Including the gradient free energy we assume the fol- 
lowing simple form for the total free energy, 

F = J dr[/ tot (0,n) + ^|V^| 2 ], (3.1) 

where / tot is given by Eq.(2.3) and C is assumed to be 
a positive constant. We suppose a planar interface at 
z = varying along the z axis and separating two phases 
a and j3. Then (f> = <f>(z) and n — n(z) are functions of z 
with </>(— oo) = 4> a and <j){oo) = <fip. In equilibrium, the 
functional derivatives of F with respect to <j) an d n arc 
constant in space. Therefore, 

h = 6F/S(f) = f'{4>) - Tgn - C<p" = constant, (3.2) 

where cf>" — d 2 (f)/dz 2 . Since u = SF/Sn = df to ta/dn is 
given by Eq.(2.7), the solute density depends on <p as 



where the coefficient A is determined by Eq.(2.9) at 
fixed n = (n). 

We multiply Eq.(3.2) by (j)' — d<j>/dz and integrate it 
over z to obtain 



Aco(<P)=co(<P)-co(<j>p) = -C(<j>'f 



(3.4) 



where uj(4>) = f—Tn—hcj) is the grand potential density in 
Eq.(2.16) satisfying w(0 Q ) = u)(<j>p) as in Eq.(2.17). From 
Eq.(2.15) and the second line of Eq.(2.16), the difference 
Abj{4>) is written as 



Aw(4>) = /(</») - f(<f> p ) - f'ififiW - 

-Tnp[e^~M -g(<f>-<Pp)-l]. 



(3.5) 



It is worth noting that the limit "f a — > can be attained 
at n = n p . From Eq.(3.5) we find 

lim Aw(<j>) = n(<j>) - [fi(0 Q )/*(^ a )]*(^), (3.6) 

where and #(0) are defined in Eqs.(2.37) and (2.38). 
From Eq.(2.45) the derivative of the left hand side of 
Eq.(3.6) with respect to <p vanishes at <f> — <f> a . There is 
no singular behavior in the interface profile even in the 
limit -f a — > 0. 

The relation Auj(cj)) = u"((f>p)((f> - 4>p) 2 /2 for = <f>p 
follows directly from Eq.(3.5). The Aw{(j>) also behaves 
as u" '{4> a ){(j) - <j> a ) 2 /2 for <j> = 4> a from Eqs.(2.15) and 
(2.17). Here u"{4>) = f"(<f>) - Tg 2 n{<j)) is the second 
derivative of Aui(4>) with respect to 4> and the correlation 
lengths £ a and £p in the bulk two phases are given by 



GC a 2 = f"{<p a )-Tg 2 n ai 
CG 2 = f"{<Pp)-Tg 2 np, 



(3.7) 



in terms of which <f>(z) — <j) a oc e~^^ a for z < — £ a and 
4>{z) — 4>p oc e~ z l^v for z > £p. The right hand sides of 
Eq.(3.7) should be positive for the stability of the bulk 
two phases. However, since f"{4> a ) ~ T/v (l — <f> a ), £ a 
becomes even shorter than a = v^/ 3 for <f> a = 1 and 
C <~ T/a. Thus the composition varies too steeply in 
phase a in the present gradient theory. However, we shall 
see (around Eq.(3.11) below) that the surface tension is 
mainly determined by the composition variation in the 
phase (3 side of an interface. 

In Alj(4>) in Eq.(3.5), the last term in the second line 
grows as exp[g(ip — <fip)] for 4> — (f)p ~> 1/g and can balance 
the other terms with increasing <f> even for very small np. 
In Fig. 7, we illustrates this behavior by plotting the sum 
of the first three terms and the minus of the last term 
in the right hand side of Eq.(3.5) separately. That is, we 
write Abj{4>) = Auj ((j)) — uj so \((j>) with 



Aoj (4>) = f{4>) - f{cj>p) 
u so \(4>) = Tnp[e 



g(<t>-<t>p) 



-f'(<Pp)(<j>-<j>p), 
gtf-fa)-!]. (3i 



n(z) = Aoexp\g<p(z)], 



(3.3) 



The Aojo((j>) is the grand potential difference without so- 
lute. As a function of 0, it increases monotonically for 
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x = 1.84 
9 = 11 

v n = 6 x 10 4 / 
A = 0.25 A> = 



£/, = 0.249 



0.2 0.4 0.6 




1 0.1 0.3 0.5 0.7 0.9 



FIG. 7: Acjo(</>) and 0J 3O i(<p) in eq.(2.59) in units of T/vo for 
X = 1.84 (left) and 2.2 (right), where g = 11, = 0.25, and 
won = 6 x 10 -4 . The solute part Wsol(</>) is negligibly small for 
4> < 0.6 here, but it grows abruptly for larger (j> and becomes 
equal to Auo(</>) at <f> = (j> a . This ensures the existence of an 
interface connecting phases a and /3. 
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FIG. 8: Left:J/(<^) in Eq.(2.56) vs (j>, which vanishes at (j> = 4> a 
and cj>i3. Right: 4>(z) vs z/a, where a = Vq and C = T/a. 
Here g = 11, cf> = 0.25, and won. = 6 x 10~ 4 . In the upper 
plates three curves correspond to (a) \ = 1-98, (b) 1.88, and 
(c) 1.84, while in the lower plates we set (d) \ = 2.2, (e) 
1.92, and (f) 1.84. In the parentheses the volume fraction -y a 
multiplied by 10 4 is written for each curve (left). 



X < 2 and exhibits positive extrema for \ > 2 (with <fip 
being outside the coexistence curve). 

In the Ginzburg-Landau scheme, the surface tension a 
is given by the z integration of the (generalized) grand 
potential density Auj(4>) + C((f>') 2 /2 including the gradi- 



ent contribution. Use of Eq.(3.4) yields 

a = / dzC^'f = y/CfJ^ f " d<t>H{J>), 

J J Aa 



(3.9) 



where we define the dimensionless function, 

H{<t>) = (2VT) 1/2 M0) " ^)] 1/2 - (3-10) 

In our theory, C is an arbitrary constant and a tx C 1 / 2 . 
In Fig.8, we plot #(0)_and <f>(z) for (a) \ = 1-84, (b) 
1.88, and (c) 1.98 with <f> — 0.46 in the upper plates and 
for (d) x = 1-84, (e) 1.92, and (f) 2.2 with 4> = 0.25 
in the lower plates. Here g = 11, vqu = 6 x 10~ 4 , and 
C = T/a. The values (a 2 a/T, n/n p ) in these cases are 
(a) (0.044,1.01), (b) (0.024,1.17), (c) (0.0021,2.25), (d) 
(0.19,6.6), (e) (0.16,7.7), and (f) (0.064, oo). In the case 
(c), the system is close to the solvent criticality and ^ 
is relatively long. In the case (f), (f> is slightly larger 
than the composition 0.249 on the oil-rich branch of the 
coexistence curve without solute, so phase separation oc- 
curs even without solute or n p = 0. The shape of H((f>) 
in Fig.8 is roughly triangular except (f), so we obtain a 
simple estimate, 



C(A<j>) 2 /^. 



(3-11) 



Since £ Q is considerably shorter than ^ , we may approx- 
imate the interface profile as 



A<f>exp(-z/£p), 



(3-12) 



in the region < z < neglecting the variation in 
the region — £ Q < z < 0. Substitution of Eq.(3.12) 
into Eq.(3.9) yields Eq,(3.11). The behaviors (3.11) and 
(3.12) hold even in the limit j a — » 0. 



B. Discontinuous appearance of a droplet due to 
surface tension: minimum droplet radius R m 

As an example, we suppose a single spherical droplet of 
phase a suspended in phase j3 in equilibrium, accounting 
for the effect of the surface tension a. The droplet radius 
R is so small that the following condition holds: 



u = ^((f> a )~/ a < 1, 



(3.13) 



where j a = 4irR 3 /3V. and $(<t> a ) ~ e^ 1 "^ . We add the 
surface free energy 47rcri? 2 to AF in Eq.(2.36) and expand 
the logarithmic factor there as ln(l + u) = u — u 2 /2. The 
resultant total free energy reads 

AF tot = V{-w la + ^Tnu 2 ) + 4naR 2 , (3.14) 
up to second order in j a . The coefficient w is defined by 



-ft(0 Q )+Tntf(0„) 



= Tv, 



\Va 



1) 

2 (x~Xp)i 



(3.15) 
(3.16) 
(3.17) 
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where n p in the second line is the minimum solute density 
in Eq.(2.46) and Xp m the third line is the minimum 
interaction parameter in Eq.(2.48). A droplet can exit 
only for w > 0, which means n > n p or x > Xp- It is 
convenient to introduce two characteristic lengths by 

R c = 2<t/w, (3.18) 
R m = (3aV/2TrTh) 1/4 y((j) a )- 1/2 , (3.19) 

in terms of which AF tot in Eq.(3.14) is rewritten as 

AFtot/47rcr = R 2 - 2R 3 /3R C + R 6 /3R^, (3.20) 

For a ~ T/a 2 we roughly estimate 

R c ~ a/(n/n p -l)~a/(x-Xp) 
Rm ~ {n p /nff\Va) 1 f 4, e- B ^-^'\ (3.21) 
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FIG. 9: (Color online) AFtot including the surface free energy 
vs R/a for (a) \ = 0-82, (b) 0.808, and (c) 0.8 from below. 
At the minimum with positive R, the droplet is stable for (a), 
marginal for (b), and metastable for (c), while n/n p = 1.080 
for (a), 1.073 for (b), and 1.067 for (c). Here R = 18.23a at 
the minimum for (b). 



If n (or x) is increased above n p (or Xp), Rm soon exceeds 
R c for large V, despite the reducing factor e -3 ' 1- ^' 4 . 

We should minimize AF to t in Eq.(3.20) as a function of 
R. For R > we require AF tot < and dAF tot /dR = 0. 
The latter condition is written as 

1 - R/R c + R 4 /R 4 m = 0, (3.22) 

under which AF tot = 47ri? 2 (2-i?/i? c )/3 < 0, so R > 2R C 
is needed. This inequality can be satisfied only for 

R c <R m /2 or w>4a/R m ,. (3.23) 

which will be evident in Fig. 18 below. If the above con- 
dition is not satisfied, there is no equilibrium droplet. As 
Rm/Rc —> 2, we have R — > R m . Here u in Eq.(3.13) 
is of order (n p /n)^ 4 [a 3 e 9< - 1 -^ /V} 1 ^ 4 and can be much 
smaller than unity for large V. Thus R m is the mini- 
mum equilibrium droplet radius in the presence of the 
surface tension. For R m /R c ^> 1, the surface free energy 
becomes negligible and 

R = {Rt/Rc) 1 ' 3 = [Ww/Ann^i^) 2 ] 1 / 3 . (3.24) 

We also note that the transition occurs for n/n p — 1 > 
a/R m from Eq.(3.16) or for (1 - 4>) 2 { X ~ Xp) > a/Rm 
from Eq.(3.17). 

In Fig. 9, we plot AFtot = AF + AnaR 2 as a function 
of R in the spherically symmetric geometry, where tf> = 
0.35, vqu = 6 x 10~ 6 , and g = 11. The cell volume 
is V = 4ttL 3 /3 £S 0.9 x 10% with L = 600a. Here 
we use Eq.(2.36) for AF and not the expansion form 
of AF in Eq.(3.14). For the three curves in Fig. 9, the 
parameter u in Eq.(3.13) is surely smaller than unity and 
Eq.(3.14) or Eq.(3.20) is a good approximation. In fact, 
at the minimum position of the curve (b) in Fig. 9, we 
have R = 18.23a, a = 0.231T/a 2 , and fi = 0.706T/v , 
while R m = 17.89a from Eq.(2.70) in this case. In the 
next section, we will present simulation results in the 
same situation for hydrophilic ions. 



C. Solvation-induced prewetting transition 

We note that Eq.(3.14) can also be used to examine the 
composition profile 4>{z) near a boundary wall at z = 0. 
Here, slightly below the precipitation curve X = Xp ( or 
n = Tin) , we may predict a first-order prewetting transi- 
tion [3(1 H3| at x = Xtr(< Xp) w ith increasing x> where 
the surface adsorption changes discontinuously. This re- 
sult is related to the experiment by Beaglehole [21| . 

Far from the wall, we assume that <j)(z) and n(z) tend 
to cj> and n, respectively. Then Eq.(3.4) holds with 

Aw(» = Q(4>) ~ Tn$(<t>), (3.25) 

where Q(tf>) and <&(</>) are defined by Eqs.(2.37) and 
(2.38). The previous Afi(^) in Eqs.(3.4) and (3.5) be- 
comes the above AQ((f>) if 4>p and np there are replaced 
by (f> and fi, respectively. Below the precipitation curve 
X = X P , Aui(4>) is positive for 4> > 4> and again approaches 
zero at cf> = 4> a as x ~ * Xp- See Fig. 7 for its behavior. 
Assuming cj)'(z) = d<fi(z)/dz < 0, we rewrite Eq.(3.4) as 

^ = -J2AM/C. (3.26) 
dz 

In addition to the bulk free energy F in Eq.(3.1), we 
assume the surface free energy, 

F s = - J dSa w (j) — ~Sa w (f) s , (3.27) 

where J dS is the surface integral at z = 0, S is the sur- 
face area, cj> s — 0(0), and a w is a parameter arising from 
the short-range interaction between the solvent and the 
surface. Here a w > if the boundary wall is hydrophilic 
and the solute is hydrophilic. The boundary condition of 
4>{z) at z = is given by 

Cty'(0) = -a,. (3.28) 
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FIG. 10: (Color online) Left: v / 2v Aui(<j))/T vs <j> with = 
0.36 and v n = 2x 10" 4 , where Acj(^) is defined in Eq.(3.25). 
For Qf w = 0.165(CT/?jo) 1 '' 2 , a first-order prewetting transition 
at a wall occurs at \ — Xtr — 1.885, where the areas (B) 
and (C) coincide. Bulk two-phase coexistence is realized for 
X > Xp = 1-8997. Right: 4>(z) for C = T/a near the wall 
before and after the prewetting transition. 



Then y / 2CAw(0) = a w at <fi — 4>s- The surface free 
energy density per unit area is then written as [3(1 [37[ 



d0[V2CAw(^) - a w ] - a w 4>. (3.29) 



Without solute, the composition increases by a w ^/C 
on the surface for \a w \ <C T/a 2 , where £ is the correlation 
length. However, in the presence of solute with j > 1, 
the adsorption on the surface can be strong with S close 
to unity for ga w ^/C > 1 even when n is very small. In 
such cases, w so \ and Awo in Eq.(2.59) are of the same 
order on the surface (see Fig. 7). In the left panel of 
Fig. 10, y/2voAu(<j))/T vs is displayed for \ = Xtr = 
1.885 and x P = 1.8997, where <f> = 0.36, n = 2xl0~ 4 , and 
a w = 0.165(CT/i;o) 1 / 2 . It demonstrates the presence of 
a first-order prewetting transition. In the right panel 
of Fig. 10, the composition profile <j>(z) is shown slightly 
before and after the transition. 



IV. HYDROPHILIC SALT 

We now treat aqueous mixtures containing a small 
amount of hydrophilic monovalent salt. We shall see that 
the charge imbalance appears only near the interface as 
an electric double layer. In the statics of the present 
problem, the role of the electrostatic interaction is thus 
to shift the surface tension slightly (see the last para- 
graph of Sec.IVA)[23]. In the dynamics, on the other 
hand, electric double layers around the droplet surfaces 
should suppress fusion of approaching droplets. Such dy- 
namical aspects are beyond the scope of this paper and 
should be studied in future. 



A. Preferential ion solvation 

We write the cation density as n\ and the anion density 
as ri2- Their total amounts are fixed as 



drni 



drri2 = Vn/2. 



(4.1) 



These densities both tend to n a /2 or np/2 in the bulk 
regions of phase a or j3 due to the charge balance. The 
electric charge density e(n\ — n-i) gives rise to the electric 
potential $ satisfying the Poisson equation, 



V • eV$ = -47re(rii - n 2 ). 

The dielectric constant e{<f) can depend on 
been assumed to be of the form, 

e{4>) = £ +£l0. 



(4.2) 
and has 

(4.3) 



Neglecting the image interaction [23J, |34| , we assume the 
total free energy F in the following form, [23| 



F = I dr 



/W0 + ^W| 2 + -^|v<rf 



i=l,2 



(4.4) 



where the third term is the electrostatic contribution and 
the terms proportional to gi represent the solvation in- 
teraction among the ions and the solvent composition. 
The Ai and A2 are the thermal de Broglie lengths of the 
two ion species. 

Let us explain the solvation terms in more detail. The 
ion chemical potentials due to solvation in aqueous mix- 
tures stro ngly depend on the ambient solvent composi- 
tion (f> [23|, [33[. We write them as nl ol (<p), where % repre- 
sents the ion species. We assume the linear form, 



A*soi( < / ) ) = Mo - T 9i& 



(4.5) 



where the first terms are irrelevant constants. This linear 
dependence is adopted to gain the physical consequences 
in the simplest manner and should not be taken too se- 
riously. In aqueous solutions, gi > for hydrophilic ions 
and gi < for hydrophobic ions. When two phases a 
and P coexist, the difference 



Am: 



/4oi(0/j) - /4 



(4.6) 



is called the Gibbs transfer free energy in electrochem- 
istry (usually measured in units of kJ/per mole), lead- 
ing to an electric potential difference across an interface, 
called the Galvani potential difference. If we assume the 
linear form (3.5), it is expressed as 



AfiU = T 9l A<t>- 



(4.7) 



For example, in water-nitrobenzene(NB) at 300K, it can 
be estimated as 17T for Na + and as 19T for CI" per 
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ion. In this mixture, the two phases are strongly seg- 
regated (Aip = 1) and we estimate g^ — 15-20 for 
small hydrophilic monovalent ions. In heavy water-tri- 
methylpyridine(3MP), the selective solvation has been 
found to induce mesoscopic phases with periodic struc- 
tures for antagonistic salt. The Gibbs transfer free en- 
ergies of hydrophilic ions for water-3MP should not be 
larger than or at least of the same order as those for 
water-NB, in view of the fact that the dielectric constants 
of 3MP and NB are 9 and 35, respectively. 

We may calculate the composition and ion profiles from 
the homogeneity of the chemical potentials h — 5F/Scf> 
and [ii = 5F/5n.i, where 



/' - cv 2 



giUi. 



(4.8) 



If Hi arc homogeneous, wc obtain the modified Poisson- 
Boltzmann relations, 

m = \np exp[gi(<f> - fo) T e($ - ®p)/T], (4.9) 



where — is for i — 1 and + is for i = 2. The potential $ 
tends to constants $ a and in the bulk regions. From 
the charge neutrality in the bulk regions, the potential 
difference between the two phases is given by 

A$ = $ Q - = T(.gi - ff2 )A0/2e. (4.10) 

Then the ion densities in the two phases satisfy 

n a /n = exp[(si + g 2 )A(j)/2}. (4.11) 

The charge density appears only near the interface for 
.9i 7^ 92 forming an electric double layer (see Fig.ll). As 
a result, the bulk phase relations in Eqs.(2.8)-(2.18) for 
neutral particles still hold for ions with the correspon- 
dence relations, 



n = n 1 +n 2 , g = (gi + ffa)/2. 



(4.12) 



For gi = g 2 we have ri\ = n 2 = n/2. 

The thickness of the electric double layer is of the order 
of the Debye screening length, X a in the phase a side and 
A/3 in the phase /? side, where 



Ak = (s K T/4im K e 



2U/2 



(4.13) 



where K = a or f3. For hydrophilic ions, A tt < holds 
from n a ^> np and A Q is longer than the interface thick- 
ness except close to the solvent criticality. We assume 
that the domain size is much longer than \p. 

We found that the surface tension of an interface of a 
precipitated domain is slightly decreased by the electric 
double layer appearing for g\ ^ g 2 . In fact, for vqu ~ 
10~ 3 , it is larger by a few % for g\ — 15 and g 2 = 7 
than for g\ = g 2 = 11, where the other parameters are 
common. 
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FIG. 11: (Color online) R/a vs \ (bold line) with hydrophilic 
ions calculated from the free energy (4.4) for (a) vqu = 3 x 
10 -4 (blue line) and (b) 6 x 10~ 4 (red line). Plotted also are 
theoretical curves from the approximate free energy (2.36) 
without the surface free energy (dotted line). Inset: Shapes 
of a spherical water-rich droplet in equilibrium for \ — 1.8, 
2.0, and 2.1 at v n = 6xlCT 4 , in the region -400 < x,y < 400 
and z — 0, where the colors represent the water composition 
according to the color bar. 
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FIG. 12: (Color online) Normalized potential e$(r)/r (with 
$/3 = 0) and water volume fraction <j)(r) (left). Normalized 
ion densities voni(r) and voU2(r) (right). Here x — 1-7 and 
eA$ = 3T. The other parameters are the same as in Fig.ll. 
In the phase /3 region near the interface, <I?(r) well agrees 
with the one-dimensional solution of the nonlinear Poisson- 
Boltzmann equation (dotted line) and there appears an elec- 
tric double layer due to nonvanishing m — ri2- 



B. Proposed experiment: a single droplet 

We performed simulation of a single water-rich droplet 
in the presence of hydrophilic ions. We have already pre- 
sented a theory for neutral solute in the same situation in 
Subsection IIF. Experimentally, wee may also suppose a 
collection of monodisperse droplets much separated from 
one another, where the volume per droplet (the inverse 
of the droplet density) should be treated as the system 
volume V. 

The simulation details are as follows. We calculated 
equilibrium profiles of <j)(r) and rij(r) by minimizing F 
in Eq.(4.3) around a spherical water-rich droplet with 
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radius R in a spherical cell with radius L = 600a. At 
4> = 0.35, we set (a) h = 6 X 10~ 4 or (b) 3 x 10~ 4 . The 
other parameters are g\ — 15, g 2 — 7, C — T/a, and 
e 2 /T = 120a. The dielectric constant depends on <fi as 
e = 40(1 + </>). Since gi > g 2 , an electric double layer is 
produced at the interface here. 

In Fig.ll, we show the equilibrium R vs x together 
with shapes of the droplet for three values of \- The 
droplet disappears due to the surface tension at R = 
18.52a and \ = 0.799 for (a) and R = 22.18a and x = 
1.656 for (b). The corresponding value of n/n p is 1.067 
for (a) and 1.102 for (b). These curves are close to those 
calculated from the approximate free energy (2.36) before 
the droplet disappearance. The characteristic features 
can also be calculated from the more approximate free 
energy (3.14) or (3.20). 

In Fig. 12, the profiles of the potential $(r) and the ion 
densities ni(r) and 712(7-) are presented for x — 1-7 in 
the case (a), which follow from Eqs.(3.2) and (3.7). Here 
4> a = 0.993 and n a = 0.352w ( ^ 1 within the droplet and 
cj)p = 0.349 and rip = 2.55 x lO^ 4 ^ 1 outside it. The 
potential $(r) relaxes with the Debye length A/3 = 11.8a 
and is well fitted to the one-dimensional solution of the 
nonlinear Poisson-Boltzmann equation (dotted line). In 
Fig. 13, the equilibrium R is displayed as a function of x 
for (a) van = 6 x 10~ 4 and (b) 3 x 10~ 4 . Remarkably, 
the droplet disappeared at R = 18.52a and x = 0.799 
for (a) and at R = 22.18a and x = 1-656 for (b) due 
to the surface tension effect. This critical radius for (a) 
nearly coincides with the droplet radius R = 18.23a at 
the minimum of the curve (b) in Fig. 9. We recognize 
that the relations of R vs x here are little affected by the 
electric double layer. Hence they are close to those of 
neutral solvent with g = 11 for the same (f> and n. 



A. Linear stability and metastability 

We first examine the linear stability of a one-phase 
state with (</>) = <\> and (n) = n. If / tot in Eq.(2.3) is 
expanded with respect to the deviations S(f> — <j> — ef> and 
Sn = n — n, the second-order free energy deviation is 
written as 

(<5/tot) 2 = \j" (0) W>) 2 - TgS<pSn + ^{Snf . (5.1) 
2 zn 

If the right hand side is nonnegative-definite, the one- 
phase state is stable or metastable. Minimization with 
respect to Sn is achieved at Sn = gn5<f>, leading to 
(<5/tot) 2 = [f"{4>) - Tng 2 ](5<f>) 2 /2. Thus the spinodal 
curve is expressed as 

f"$)-Tng 2 =0. (5.2) 

See the left panels in Fig.l for this spinodal curve. One- 
phase states are linearly stable outside this curve. 

In the presence of the Coulomb interaction, the spin- 
odal curve is still given by Eq.(5.2) with the replacement 
g = (gi + g 2 )/2 in Eq.(4.12) under the condition 7 P < 1. 
Here 7 P is a parameter representing the degree of solva- 
tion asymmetry between the cations and anions (which 
should not be confused with the volume fraction 7 Q of 
phase a). For monovalent ions it is given by 

7p = (T/167rC) 1 / 2 | ffl - 52 |. (5.3) 

For most hydrophilic ion pairs, we should have j p < 1. 
On the other hand, the reverse condition 7 P > 1 can be 
realized for antagonistic salt composed of hydrophilic and 
hydrophobic ion pairs [H, [H, 23] . Addition of such a salt 
leads to a decrease in the surface tension and mesophases 
formation. In this paper, we consider hydrophilic ions 
with 7 P < 1 for simplicity. 



V. HOMOGENEOUS NUCLEATION FROM 
ONE-PHASE STATES 

We start with homogeneous one-phase states without 
impurities (other than the solute or the salt under con- 
sideration). We assume that the boundary walls are hy- 
drophobic and formation of a water-rich wetting layer is 
suppressed. In this situation, the precipitation predicted 
in this paper can be realized via homogeneous nucleation 
for n > n p or x > Xp- Therefore, we calculate the nucle- 
ation rate in the initial stage of nucleation. It is not much 
affected by the electrostatic interaction for hydrophilic 
ions and we treat neutral solute in the following. We 
also neglect temperature inhomogeneity and treat x as a 
homogeneous constant. 

It is worth noting that water droplets can easily be 
produced around hydrophilic ions in metastable gas mix- 
tures containing water vapor. Here ions play the role 
of nucleation seeds on which hydration-induced conden- 
sation is favored. A Ginzburg-Landau approach to this 
problem was also presented [351 ] . 



B. Droplet free energy and nucleation rate 

In the early stage of nucleation, the volume fraction 
7 Q of the new phase is very small and the droplets may 
be treated independently. Hence we consider a single 
droplet with radius R in the classical nucleation theory 
[32| . The composition and the solute density far from 
the droplet are <f> and n. A key quantity here is the 
free energy AF(R) needed to create a droplet of the new 
phase. 

To derive AF(R), we use the total free energy density 
/tat in Eq.(2.3). It is equal to / tot = ftot(<t>,n) in 
the initial homogeneous state and remains so far from 
the droplet, while it is equal to / t " t = /tot (0a, n a ) inside 
the droplet. Here <f> a and n a are the composition and 
the solute density inside the droplet, respectively. The 
change of the total free energy is 

AF(R) = -jR 3 [f t a ot - /tot] 

+ / dr[f tot (4>, n) - / tot ] + AiraR 2 . (5.4) 

Jr>R 
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In the right hand side, the first term is the contribution 
from the droplet interior, the second term is that of the 
droplet exterior, and the last term is the surface tension 
term. Outside the droplet we assume small deviations of 
4> and n from the initial values <f> and n so that we set 
/tot(0, n) - /tot = h{<t> - 4>) + fi(n - n) , where h and p, are 
the values of h in Eq.(2.6) and /x in Eq.(2.7) in the initial 
metastable state. Further we note that the integrals of 
(j> — </> and n — n in the whole space vanish from the 
conservation relations (2.5), so that their space integrals 
outside the droplet are equal to — (47r/3)i? 3 (0 Q — 4>) and 
— (47r/3)i? 3 (n a — n), respectively. Therefore, it follows 
the standard form [32j |. 

AF(R) 



Air 



-R 3 w + 4itR 2 <t. 



(5.5) 



The coefficient w represents the degree of mctastability 
and is of the form, 



W = /tot - /tot - h(4> a ~ 

= n(<p a )-Tm(<p a ), 



H(n a - n) 



(5.6) 



where Sl(4> a ) and &((f) a ) are given in Eqs.(2.37) and 
(2.38). From the second line of Eq.(5.6), this w coin- 
cides with w introduced in Eq.(3.15). It is proportional 
to n — n p and x ~ Xp as m Eqs.(3.16) and (3.17), where 
n p and x P are defined in Eqs.(2.46) and (2.48). 

The critical radius is R c — 2a /w as in Eq.(3.18). 
R = R c , AF(R) attains a maximum given by 



At 



(5.7) 



In the classical nucleation theory, critical droplets appear 
with a constant rate I per unit time and per unit volume. 
It is called the nucleation rate and is of the form, 



I = I exp(-F D /T), 



(5.8) 



where the coefficient Iq depends on the dynamics and 
will be estimated in Eq.(5.22). For weal metastability, 
the nucleation barrier F c /T grows as w~ 2 and / is very 
sensitive to w. From Eqs.(5.7) and (5.8) the nucleation 
barrier is expressed as 



F c /T = A c /(n/n p - l) 2 
= B c /( X - Xp ) 2 . 

The two coefficients A c and B c arc defined by 

A c = 167rcr 3 /3ft 2 , 

B c = 16^ 3 ^/3T 2 (cA tt - 0) 4 , 



(5.9) 
(5.10) 



(5.11) 
(5.12) 



where A c depends on x and 4> and B c on n and <f>. We 
may use Eq.(5.9) or Eq.(5.10) depending on whether n 
or X — x{T) is varied. 

In Fig. 13, we plot R c vs n/n p — 1 at fixed x m the left 
and R c vs \ ~ Xp at fixed n in the right. It much exceeds 
the molecular size a for weak metastability, where n/n p — 
1 or x — Xp is small. In Fig. 14, we plot A c at fixed x 
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FIG. 13: (Color online) Critical radius R c in Eq.(3.18) as a 
function ofn/n p -lfor ( X ,(j>) = (1.8, 0.1), (1.2, 0.1), (1.2, 0.3), 
and (1.8, 0.2) from above (left) and as a function of x ~ Xp 
for (10 4 v n,(j>) = (2.0.17), (0.7,0.1), (2,0.22), and (0.7,0.15) 
from above (right). Here g = 11 and C — 2T/a. 
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FIG. 14: (Color online) A c at fixed % in Eq.(5.11) (left) and 
B c at fixed n in Eq.(5.12) (right) as functions of <j> for g = 11 
and C = 2T/a. For \ > 2 (left), there is no barrier or A c — !> oo 
as 4> approaches the equilibrium coexistence composition. 



and B c at fixed n as functions of (j) for g = 11. These 
coefficients are considerably larger than unity for 4> ^ 0.2, 
but become small with increasing </>. In these figures we 
set C — 2T/a. Here notice the relations R c oc C 1 / 2 , 
a cx C 1 / 2 , and F c oc C 3 / 2 . The coefficients A c and B c in 
Eqs.(5.11) and (5.12) depend on C as cr 3 oc £ 3 cx C 3 / 2 . 
For aqueous fluids with the hydrogen bonding network, 
a larger value of C might be more appropriate than in 
Figs. 13 and 14. See the summary for more discussions 
on the choice of C. 

Nucleation experiments have been performed precisel 
on near-critical binary mixtures (without salt) |32l |3j 
libT ]. where appreciable droplets become observable for 
Fc/T < 50. From Fig.14, we find F c /T < 50 for 
n/n p > 1.2 at ip — 0.1 In our case, once n > n p or 
X > Xpj the nucleation rate should soon becomes large 
enough for usual observations of droplets with R > R c . 



C. Weak metastability and droplet growth 

We may derive an approximate expression for w in 
terms of <p and n for weak metastablilty. To this end, we 
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consider an equilibrium reference state, where two bulk 
phases with (4>,n) — {4> a -,n a ) and (^g,^) are separated 
by a planar interface. We assume that the composition 
and the solute density (0 a ,n a ) in the reference phase a 
are the same as those within the droplet. However, those 
(cj)p,rif)) in the reference phase (3 are slightly different 
from those (</>, n) in the initial metastable state. The 
initial homogeneous deviations are written as 

8<f> = <}> — (j>fj, 8n = n — np. (5.13) 

From Eqs.(2.16) and (2.17) we obtain 

f t a ot - /f ot - h CK Acf) - McxAn = 0, (5.14) 

where f£ t = f to t((f>K ,n K ) (with K = a,/3)). The h cx 
and the values of h in Eq.(2.6) and \x in Eq.(2.7) 

in the reference state (which are simply written as h and 
/i in Sec. II). Further we note the relation, 

/tot - /fot - h cx Scf> - n cx 5n = 0, (5.15) 

which is valid to first order in 5(f> and Sn. From the first 
line of Eq.(5.6) we may eliminate / t ^ t and /tot to obtain 
the desired expression, 

w = (h - h cx )A<f> + (fi - (U cx )An, (5.16) 

where we have set <p a — <j) = A(f> and n a — h = An. From 
Eqs.(2.6) and (2.7) the initial deviations of the chemical 
potentials are expanded as 

h-h cx = fp54> - Tg8n 7 

fi- /J, cx = TSn/n/3 -TgSifi, (5.17) 

where /^' is equal to d 2 f/dcf) 2 at cf> — <f>p. Therefore, 

f = (#^r " .9An)50 + (e» A * - 1 - gA<f)5n. (5.18) 

In Appendix B, we will derive the dynamic equation 
for the droplet radius R(t) for weak metastability. It is 
of the standard form, 

-m=R{ w -R)> (5 - 19) 

where A is a kinetic coefficient defined in Eq.(BlO) be- 
low. For dilute solute, it may be related to the mutual 
diffusion constant of the mixture D m ((f>) by 

A = D„S)/f"m^) 2 - (5.20) 

For R = R c we have dR/dt = T C (R - R c ), where 

T c = 2aA/R 3 c . (5.21) 

The inverse Tj 1 is the time scale of near-critical droplets. 
The coefficient Iq in the nucleation rate (5.8) is deter- 
mined by the droplet dynamics at R = R c (with appro- 
priate thermal noises added) and is estimated as [32| 

h - rye 3 , (5.22) 



where £ is the correlation length in the initial metastable 
state. 

For g = our results tend to those for incompressible 
binary mixtures without salt [HI, HH, H3- There, we use 
Eq.(5.22) and the relation w - 2a/ R = /^'(A</>) 2 [A S - 

2do/ R] in terms of the supersaturation A s ~ 5<j)/ A0 and 
the capillary length do = 7 / fp(Acj)) 2 . 



VI. HETEROGENEOUS NUCLEATION ON 
HYDROPHILIC COLLOID SURFACES 

In this section, we examine adsorption, wetting, and 
precipitation on a colloid surface. Here the colloid sur- 
faces are hydrophilic, while the boundary walls are hy- 
drophobic as in the previous section. The colloid density 
n co is so small that each colloid particle may be treated 
independently and the composition profile around it de- 
pends only on the distance r from its center. The effective 
cell volume for each particle is V = n~^. For simplicity, 
we assume that the colloids are neutral without surface 
charge and the correlation length £ outside the colloids 
is shorter than the colloid radius d. 

In Sec.IIIC, we have already found a prewetting tran- 
sition on a planar wall for x slightly below \ p . Further- 
more, we shall see that a wetting layer much thickens 
slightly above Xp due to precipitation. 



A. Simulation results: discontinuous and 
continuous transitions 

First, we present our numerical results for a hy- 
drophilic particle with radius d — 15a or 25a placed at the 
center of a spherical cell with radius L = 10 3 a. The cell 
volume is V = 47rL 3 /3 and the colloid volume fraction is 
co i = (d/L) 3 , so <^coi = 3.4 x 1(T 6 and 1.6 x 1CT 5 for 
d = 15a and 25a, respectively. We suppose hydrophilic 
ions using the model in Sec. IV, though the electrostatic 
interaction is not essential here. The other parameters 
used are </> = 0.36, v n = 2 x 10~ 4 , g\ = 9, g 2 = 13, 
C = T/a, e 2 /T = 120a, and e x = e a = 40. The correla- 
tion length £ far from the surface is shorter than d. In 
fact, it is 1.39a at x — 1-9- 

For the composition <f> and the ion densities «i and «2, 
the total free energy is Ftot = F+F s , where F is given by 
Eq.(4.3) and F s is the surface free energy in Eq.(3.27). 
The composition <j)(r) obeys Eq.(3.2) at a constant h. 
The boundary conditions at r — d and L are given by 

C(j)'{d) = -a w , <t>'{L) = 0, (6.1) 

where 4>'(r) = d<p(r)/dr. The wetting parameter a w is 
set equal to 0.165T/a 2 as in Fig.10. 

In Fig. 15, we plot the preferential adsorption given by 

r= f dr[<f>(r) - <f>(L)}, (6.2) 
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FIG. 15: Preferential adsorption V on a colloid surface in 
Eq.(6.2) as a function of y on a semi-logarithmic scale at <f> = 
0.36. For d = 25a (upper curve), after a prewetting transition 
at x — 1-873, the crossover to precipitation is continuous. For 
d — 15a (lower curve), after a prewetting transion at x = 
1.898, the transition to precipitation is discontinuous at y — 
1.918. At the points (a)-(h), the corresponding composition 
profiles are given in Fig. 16. 



for the two diameters d = 25a and 15a. The prewet- 
ting transition discussed in Sec.IIIC occurs at x = 1-873 
for d = 25a and 1.893 for d = 15a. Remarkably, for 
X > Xp — 1-8997, T increases continuously for d = 25a, 
but discontinuously for d = 15a. In Fig. 16, we show the 
composition 4>{r) for four x across the transitions and 
around the crossover. It changes from a thin to thick 
wetting layer continuously (but abruptly) for d = 25a 
and discontinuously for d = 15a. In Fig. 17, we show the 
free energy change AF to t = -Ftot — -Ftot as a function of 
X, where F® ot is the total free energy for the homoge- 
neous state (realized for a w = 0). Around the first-order 
transitions there appear two branches, where Eq.(3.2) is 
satisfied, and on the equilibrium branch F tot takes a lower 
value. In particular, for d = 15a, two branches of adsorp- 
tion and precipitation appear around x — 1-918 > Xp- 

The prewetting transition point below Xp depends on 
d as in Fig. 15. Furthermore, it is slightly affected by the 
electric double layer for g\ ^ 71%. In fact, for g\ = gi = 
11, it disapperas and the prewetting transition occurs at 
X = 1-907 for d = 25a and at X = 1-912 for d = 15a, 
while it was at x = 1-885 for a planar wall in Fig. 10. 



B. Theory of wetting on a colloid surface 

Before the layer thickening due to precipitation and for 
d ^> £, A.Ftot is due to the adsorption expressed in terms 
of J" ad in Eq.(3.29) as 



AFtot =4^d 2 .F ad (<M), 



(6.3) 



We explicitly write the dependence of J-ad on the com- 
position <f> and the solute density n far from the surface. 
In Fig. 17, J-ad only weaky depends on x- After precipi- 
tation in the range x > Xpi a thick wetting layer appears 
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FIG. 16: (Color online) Composition <j)(r) vs r/a outside a 
colloid r > d. For d = 25a (left), X is (a) 1.87, (b) 1.875, (c) 
1.915, and (d) 1.92. For d = 15a (right), X is (e) 1.896, (f) 
1.899, (g) 1.917, and (h) 1.919. 
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FIG. 17: (Color online) Free energy change AFtot due to 
water accumulation around a colloid as a function of x f° r 4> = 
0.36. For d = 25a (left), after the prewetting transion between 
(a) and (b), the crossover to precipitation is continuous but 
there is an abrupt slope change at (c). For d = 15a (right) 
after the prewetting transion between (e) and (f), there is a 
discontinuous change between (g) and (h). Shown also is the 
approximate free energy F we t in Eq.(6.6) (dotted line). 



in the region d < r < R and A_F t ot is expressed as 

AF tot = 47rd 2 F ad (0 Q ,7i Q ) + F wot . (6.4) 

The first adsorption term is determined by the compo- 
sition <p a and the solute density n a in the surrounding 
thick layer, while the second one is the contribution in 
the thick layer. Figure 17 indicates that F we t decreases 
dramatically with increasing x- 

To construct a simple theory of F wctl we assume the 
condition ^(4> a )"f a <C 1 in Eq.(3.13). We treat AF tot in 
Eq.(2.65) as F WQt . In the present situation, the volume 
fraction of phase a is given by 

47r(i? 3 - d 3 )/3V = (R 3 - d 3 )/L 3 . 



la 



(6.5) 

Here w in Eq.(3.15) is positive. In terms of the two 
lengths R c oc w^ 1 in Eq.(3.18) and R m in Eq.(3.19), we 
may express -F we t as 

1 -(R 3 -d 3 ) 2 



^ = R 2 -d 2 -—(R 3 -d 3 ) 
4 tt(t 3R C y ' 



3i?i 



— R„ 



D 3 )2 /3_ D 2_^ g+ l g 2 



(6.6) 
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The first line follows from Eq.(3.20) if R 2 and R 3 there 
are replaced by R 2 — d 2 and R 3 — d 3 , respectively. Thus, 
as d — > 0, this F wet tends to AF to t in Eq.(3.20). In the 
second line we introduce the order parameter, 



(R 3 - d 3 )/R. 



3 

m ! 



(6.7) 



where D and E are dimensionless parameters defined as 

D = d/R m , E = R m /R c = wR m /2a. (6.8) 

We minimize F wct as a function of R or q. In Fig. 17, this 
-Fwet is plotted (in dotted lines) and is compared with 
AFtot- In Fig. 18, we display q in the D-E plane, which 
shows how a thick wetting layer appears. For D = 
(or d = 0) we have the curve of R 3 /R 3 1 vs R m /R c for a 
droplet with radius R (see Sec. IFF for its theory). 

For d — R » £, the adsorption free energy 47rd 2 J 7 a d 
is small as compared to F wot . Then we may examine 
the formation of a thick layer by minimizing Fwet- For 
q -C D 3 , Fret is expanded in powers of q as 



AnaRl, 



-(- 

3 y D 



E)q+('S 



D 4> 9 



4<f 
81L> 7 



(6.9) 



By setting the coefficient of q 2 equal to zero, we find a 
tricritical value of D given by 



Dtri = 3~ 1/4 = 0.760. 



(6.10) 



For D = d/Rm < Ari, the second term (oc q 2 ) in the 
expansion (6.9) is negative and the transition takes place 
discontinuously on the transition line E — E tI (D), where 
F tr -> D~l as D — > Ari and F tr — > 2 as D — > 0. On the 
other hand, for D = d/R m > Arb a thick wetting layer 
appears continuously as a second-order phase transition. 
The parameter q is nonvanishing for 



E> D~ 



or d> R c 



.11) 



where q oc (F — D 1 ) 1 ' 2 . Thus a tricritical point is at 

(D,E) = (Ari, Art) OTat K w ) = (Ari#m,2a/Ari#m) 

for precipitation on a colloid surface. In Figs. 15-17, we 
have i? m = 26a, so that D — 0.96 > Ari for d = 25a and 
D = 0.58 < Ari for d = 15a. Our numerical results are 
consistent with the prediction from Eq.(6.6). 



VII. SUMMARY AND REMARKS 

In summary, we have examined solvation-induced pre- 
cipitation in aqueous mixtures with hydrophilic or hy- 
drophobic solute in the limit of large g, which repre- 
sents the composition-dependence of the solvation (the 
strength of the preferential solvation). 

In Sec. II, we have developed a thermodynamic theory 
of two-phase coexistence in mixture solvents with solute, 
numerical analysis of the two-phase coexistence, a the- 
ory of precipitation in the asymptotic limit j > 1, and 




FIG. 18: (Color online) Precipitated layer volume 4-7r(i? 3 — 
d 3 )/3 divided by 47ri?f n /3 in the plane of d/Rm and Rm/R c - 
It is calculated from the approximate wetting free energy F wot 
in Eq.(6.6) without the adsorption free energy 47rd 2 J r ad . The 
transition is discontinuous for d/R m < 0.760 on the perpen- 
dicular surface (in red). 



a theory in the dilute limit of precipitated domains. Re- 
markably, the precipitation curve x = Xp extends far 
below the coexistence curve in the X'4 1 plane. In Sec. Ill, 
introducing the gradient free energy, we have calculated 
the surface tension er, performed stability analysis of a 
precipitated droplet including the surface tension, and 
predicted a prewetting transition below the precipitation 
curve. In Sec. IV, we have presented a theory for hy- 
drophilic ions including the solvation and electrostatic 
interactions. Finally, we have investigated homogeneous 
nucleation in Sec.V and heterogeneous nucleation on a 
collod surface in Sec. VI. We have also derived the evo- 
lution equation for the droplet radius R(t) from the dy- 
namic equations for <j> and n in metastable states in Ap- 
pendix B. 

We have presented a number of predictions. They are 

(i) droplet appearance due to preferential solvation as 
in Figs. 1-3, which has a minimim size i? m in Eq.(3.19), 

(ii) droplet size change with varying x or T as in Fig. 11, 

(iii) prewetting transition slightly below the precipita- 
tion curve, (iv) nucleation slightly above the precipitation 
curve and far outside the solvent coexistence curve, and 
(v) heterogeneous nucleation on a collod surface, where 
two-step first-order wetting transitions can occur as in 
Figs. 15-18. To confirm these effects, systematic experi- 
ments are needed, where the temperature, the water vol- 
ume fraction, and the salt amount are controlled. Partic- 
ularly, dynamic light scattering is informative to detect 
emergence of droplets and wetting on colloid surfaces. 
Wetted colloids should move with a smaller diffusion con- 
stant oc R~ x . 

In our theory, the molecular volumes of the two com- 
ponents are assumed to be given by the common vq in 
Eq.(2.4). However, they can be very different. For ex- 
ample, those of D2O and 3MP (the inverse densities of 
the pure components) are 28 and 168 A 3 , respectively. 
Moreover, the coefficient C in Eq.(3.1) of the gradient 
free energy remains an arbitrary constant and the nucle- 
ation rate calculated in Sec.V is proportional to C 3 / 2 . 
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Therefore, our theory is still very qualitative. For pure 
water, the observed surface tension outside the critical 
region 1 — T/T c > 0.1 fairly agrees with the calculated 
surface tension from the van der Waals model with the 
gradient free energy density 10Ta 5 |Vn| 2 /2 [IHHlJ], where 
n is the density and a = 3.lA is the van der Waals radius. 

In future, we will investigate the wetting transition 
on charged walls, rods, and colloids and the solvation- 
induced colloid interaction, which can be much influenced 
by the ion-induced precipitation. The degree of ioniza- 
tion should also be treated as a fluctuating variable sen- 
sitively depending on the local environment (29j . 
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Appendix A: Two-phase coexistence for \ > 2 
at very small n 

For x > 2 two-phase coexistence is possible without 
solute. Let </>° and (jPp be the equilibrium volume frac- 
tions in the two phases without solute. Here we calculate 
the deviations of the volume fractions 5(f) a = 4> a — 4> a and 
5<f>[} — <pp — <f>g to linear order in the solute density n in 
the mean field theory. We show that the linear theory 
holds only for very small n. 

The n a and np are still given by Eqs.(2.8) and (2.9). 
From Eqs.(2.15) and (2.17) the chemical potential differ- 
ence h is expanded as h = — TAn/A(f> so that 

6cb a = Tn a (gA(f> - 1 + e' 9 ^)/ f" A(j>, 

--Tn e {gA^ + l-e 9 ^)/f"A^ (Al) 



where /" = f"(<fi%) = /"(0a) for the symmetric free 
energy density in Eq.(2.4). Here /"(</>) = d 2 f/dcj) 2 . 

As the solvent critical point (\ c = 2) is approached, 
the inequality Acf> < 1/j holds eventually to give 



Tng 2 A<t>/2f". 



(A2) 



If we set e = \ — 2(oc 1 — T/T c ), the above deviation is 
of order v^ne^ 1 / 2 . The linear approximation holds when 
Ha < 0° - 1/2 = A<f)/2. Thus the upper bound of the 
linear regime of the solute doping is very small as 



n < e 3/2 /g 2 v . 



(A3) 



In the left panels of Fig.l, cj) a — <f>p is not small even for 
X = 2 due to the nonlinear solute effect. 

Away from the critical point, the relation gA<p 3> 1 
holds for large g. In this case, we find 



so that S(j) a = -(gA(/))6(p/3 > 0. From Eq.(2.13) <5</> Q (oc 
n a ) increases steeply as (f> is decreased to </>^ (or as 7 Q is 
decreased). This tendency can be seen in the left panels 
of Fig.l. The maximum of n a is he 9A ^ . The condition 
that 6(j) a is much smaller than unity is written as 

n< f"e- 9 ^/Tg (A5) 

If Eq.(A3) or Eq.(A5) does not hold, cj) a nonlinearly de- 
viates from <jP a with respect to n. 

Appendix B: Droplet growth 

For weak metastability, we derive the droplet-evolution 
equation (5.19) in the neutral solute case. The simplest 
dynamic equations for the composition <p(r, t) and the 
solute density n(r, t) are given by the diffusive equations, 



at 
dn 
~dt 



LV 2 (h-CV 2 (/>), 
D S V • nV ■ (jt/T), 



(Bl) 
(B2) 



where L is the kinetic coefficient for the composition, D s 
is the solute diffusion constant, and h and are given 
in Eqs.(2.6) and (2.7). Around a spherical droplet with 
radius i?, all the quantities depend on the distance r from 
the droplet center and the time t. 

Slightly outside the droplet surface r — R > the 
gradient term in h is negligible, where £ is the interface 
thickness. When the droplet growth or shrinkage is slow, 
we may use the quasi-static approximation: 



h = h + (h R - h)R/r, 
/i = fi + (i-iR - p.)R/r, 



(B3) 



where h and /i tend to h and jl far from the droplet and 
to hji and /ir near the droplet surface. The conservations 
of 4> and n at the interface yield 



9R AI L n 7N 

~~dt = ~R ~ 

OR. _D 

—An = n-(/x R - M ) I 



(B4) 



where A<f> = <f> a — cf>p and An = n a — np as in Eq.(2.13). 

In the interface region \r — R\ < £, the generalized 
chemical potential //— CV 2 (f) including the gradient term 
may be treated as a constant Hr, so that 



or z r or 



(B5) 



The solute chemical potential /i is also a constant hr 
near the interface. We set <j> — <fii n t(r — R) + 5<j> and 
ti = rii n t(r — R) + Sn in the interface region. Here 4>i n t(z) 
is the one-dimensional interface solution of Eq.(3.4) and 
nint{z) — n ( 3e ff< ^ int ( z - l_ *' 3 ). These represent the reference 
profiles and the corresponding h and ijl are written as h cx 
and fj, cx as in Eq.(5.14). To linear order in 8(j) and 5n, we 
obtain (1r — fi cx = T5n/n ln t — Tg8<j) or 



Tgnjf", Sfo Si -Tn a /f"A(t>, 



(A4) 



Sn = n int [(fi R - ^ CX )/T + g5(p], 



(B6) 
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in the region \r — R\ < £. We may linearize Eq.(B5) as 

d 2 



f" 



Tg 2 n int - C 



dr 2 



R 



C 0int ~ 9(^R - McxHnt, 



(B7) 



Then we multiply Eq.(B7) by (j>[ nt = d<p- m t (r — R)/dr and 
integrate over the region \r — R\ < £ to obtain 



2a 



(Mfl. - Mcx)An, 



(B8) 



where the first term in the right hand side of Eq.(B7) 



does not contribute and use has been made of Eq.(3.9) 
and gni n t<?!>[ nt = drii n t/dr. We rewrite Eq.(B8) as 

(h R - h)Acj) + (fi R - fi)An ^^-w, (B9) 

where w is given in Eq.(5.16) in the linear form. Using 
Eq.(B4) we eliminate hn — h and hr — P< from Eq.(B9) 
to obtain the droplet evolution equation (5.19) with 

1/A = {A4>) 2 /L + (An) 2 /D s n. (BIO) 

where the second term in the right hand side is negligible 
for dilute solute, leading to A = L/(A4>) 2 and Eq.(5.20). 
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